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Starting from the entropy production being invariant under time reversal, one can (i) easily 
proof, and understand, many aspects of the linear Onsager relations and (ii) deduce the result that 
all quadratic Onsager coefficients for hydrodynamic fluxes vanish. 
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Despite its usefulness, simplicity, and thorough experi- 
mental verification, the Onsager symmetry relation poses 
a number of vexing questions - with respect to the ambi- 
guity of its prescription, to its derivation and to its non- 
linear generalization. Consider the temporal evolution 
of the entropy density, s — W ■ f — TZ/T, with T denot- 
ing the temperature and f the entropy current; TZ is the 
source term, or entropy production, usually given as a 
sum of products: TZ = biBi. Expanding the fluxes 
bi in the forces Bi to linear order, bi = LijBj, the 
Onsager relation relates the kinetic coefficients pair- 
wise, Lij = tip{BiBj) Lji, where tip(^) stands for the 
time inversion parity of the quantity A. The problems 
with this simple sounding prescription are: 

(1) On the most elementary level, given a product 
biBi, it is not always clear which of the two is the force - 
the time inversion property of which is so important for 
the Onsager relation. This has been the main thrust of 
an eloquent and polemic criticism ||]. 

(2) A step further, the expression TZ — J2i ^i^i is by no 
means unique, and the question naturally arises whether 
the Onsager relation still holds for TZ = J^i biBi, where 
bi and Bi are linear combinations of bi and Bi. 

(3) In the original proof by Onsager, the fluxes are 
given by the temporal derivative of the state variables. 
This is not the case with hydrodynamic fluxes such as 
the entropy current f above; and it is not clear whether 
the symmetry relation holds here without further quali- 
fications. On the other hand, we have no hint whatever 
of any invalidation experimentally; and there have been 
successful efforts to bridge this gap for specific exam- 
ples But a generally valid, clean-cut criterion for 
deciding when a given bi is a legitimate flux has proved 
elusive so far. 

(4) The validity of the Onsager relation is explicitly 
limited to linear terms, bi = LijBj with con- 
stant. Nonlinear generalization appears difficult, and has 
focused on the question whether the symmetry relation 
such as it is holds nonlinear ly, if Lij depends on Bj ^ . 

In addition to these problems of the prescription, the 
Onsager relation has been the one result of the thermo- 
and hydrodynamic theory that could not be derived and 
justified within the framework of macroscopic physics. 
Moreover, the knowledge required for the excursion to 



microscopic dynamics and statistical mechanics is so ad- 
vanced that many textbooks prefer to only state the re- 
sults; others clad it in formidable notations and tuck it in 
the last few chapters of a thick book; the more convincing 
ones still fight the appearance of being at odds with the 
ultimate simplicity and universal validity of the Onsager 
relations. As a result, the knowledge about this relation 
has remained incommensurate with its fundamental sig- 
nificance and sheer usefulness; even the general awareness 
of it has suffered. 

So it must be a welcome discovery that there is in 
fact a simple and macroscopic derivation for the On- 
sager relation. Starting from the postulated invariance 
of the entropy production TZ under time inversion, only 
a few simple, algebraic steps are needed to conclude the 
proof, which includes the hydrodynamic fluxes and is eas- 
ily extendable to nonlinear terms. It also eliminates the 
above ambiguities, showing that a sufficient condition for 
a bona fide force is an unambiguous tip; and it specifies 
the most general linear transformation that leaves the 
Onsager relation undamaged. One such transformation 
actually reduces Lij to the unity matrix. 

Depending on sentiments, some will not want to call 
this a proof, as a new postulate is introduced. Neverthe- 
less, since the invariance of TZ concisely summarizes all 
the above points, it is - if true - certainly a much bet- 
ter starting point when thinking about the Onsager rela- 
tions, and when presenting them in books or lectures, as 
it should greatly facilitate the understanding and teach- 
ing of this chapter of irreversible thermodynamics. We 
should therefore at least accept it as a preliminary yet 
powerful hypothesis, while working to provide a proof, 
more rigorous than the considerations below - or else, 
show where it may conceivably fail. 

In the following, the case for the high plausibility of 
TZ's invariance, tip(7^) = +1 is presented first. It is then 
employed, without further input, • to derive the Onsager 
relation, • to obtain the most general linear transforma- 
tion, • to include hydrodynamic fluxes, and • to consider 
nonlinear ramifications, showing especially that all sec- 
ond order Onsager coefficients for hydrodynamic fluxes 
vanish. This is good news, as hydrodynamic theories — 
such as the Navier-Stokes equations, the superfluid hy- 
drodynamics of '^He [||, the nematodynamics 0, and the 
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hydrodynamic Maxwell equations |g| — are all an order 
more accurate then hitherto thought. 

Imagine particles in a box occupying only one half of 
the volume. Their velocity distribution is Maxwellian, 
with zero mean. Given these initial conditions, the parti- 
cles will spread out quickly, and produce entropy on their 
way. If the time runs backward, x ^ v —>■ —v for each 
particle, and no macroscopic change in the initial con- 
ditions takes place. By pure phase space arguments we 
again expect the particles to spread. The rate of spread- 
ing, and with it the rate TZ of entropy production, must 
remain unchanged if the microscopic dynamics remains 
unchanged, if it is reversible. So, micro-reversibility leads 
straight to the invariance of TZ. 

This simple scenario is of course a special case, but 
the conclusion should be generic, as the crucial element 
here, the reversibility of the dynamics and initial condi- 
tions, holds widely: The first is believed to be rigorously 
true, the second presupposes local equilibrium, (or par- 
tial equilibrium where only a few slow macroscopic de- 
grees of freedom are not in equilibrium) . Off local equilib- 
rium, say if we had more particles moving inwards than 
outwards initially, there will be a transient, ballistic pe- 
riod in which the cloud of particles will shrink, with the 
expansion resuming only afterwards. Reversing the time, 
the cloud will expand immediately. However, this tem- 
poral asymmetry disappears on a coarser time scale, on 
which the initial period is infinitesimally short. And only 
on this scale should the spreading be related to the (in- 
trinsically coarse-grained) entropy production TZ — which 
is zero in the transient range (no loss of information), 
and retains its invariance. 

If there is an external magnetic field H, micro- 
reversibility is restored if the field is flipped with the time 
reversal operation. The same is true with other macro- 
scopic quantities that are odd under time inversion. 

Summarizing, we conclude that starting from any im- 
probable, nonequilibrium state, the entropy increases for 
statistical reasons, in either directions of time. Given 
micro-reversibility, the rate of increase is the same. The 
first sentence implies TZ > 0, the second tip(7?.) — +1. 

To understand clearly the implication of these two ex- 
pressions, consider s — TZ/T and its solution s{t). Choose 
a point to in time for a system that is off-equilibrium, 
7^ 7^ 0. Note tip(s) = -1 and tip(r) = +1. If 
tvp{TZ) = —1 were true, the solution s{t) would be un- 
changed under the inversion of time, though the same 
curve is now traced backwards. With tip(7?.) ~ +1, in- 
verting the time at to reflects s{t) at to. Retaining only 
the portion t > t^ of s{t), the two portions for either di- 
rections of time meet at to and form a cusp there. With 
7^ > 0, the cusp points downwards. Accepting 7^ > 
yet allowing TZ — "^^biBi — LijBiBj to contain 
terms for which tip{BiBj) = —1, results in an asymmet- 
ric cusp, still pointing downwards but with the left arm 
(depicting the time reversed solution) being flatter. 

Two final remarks. First, the proof below is such that 



it may be traced backwards, starting from the validity of 
the Onsager relation to demonstrate the invariance of TZ. 
So, at least for the linear regime, close to global equilib- 
rium, this invariance may be considered beyond dispute. 
At the same time, none of the considerations above makes 
any reference to how removed from global equilibrium the 
initial conditions may be. There is therefore no reason 
why TZ should not remain invariant including nonlinear 
Onsager coefficients. The half filled space is certainly far 
off equilibrium, not accountable by a linear expansion in 
density gradients. Yet this is our most straight forward 
example of an invariant TZ. 

Second, when accounting for dilute systems with the 
kinetic theory, higher order gradient terms are inherently 
related to higher order temporal derivatives. The latter 
account for inertial effects (or the dynamics of additional 
variables) , and lead to "temporal slips" such as the asym- 
metric spreading discussed above. In hydrodynamic the- 
ories of dense systems, these two expansions are quite 
independent. And an inclusion of higher order Onsager 
coefficients does not drive the system off local equilib- 
rium. A comparison between the results of both theories 
must therefore be executed with circumspection 

We proceed to consider the ramifications of an invari- 
ant TZ, first in the context of non-hydrodynamic, global 
and spatially uniform variables. With Ui a complete set 
of such variables, and the entropy S{ai) maximal in equi- 
librium, we have (i) Ai = dS/dai — 0, and (ii) hi = 0. 
Slightly off equilibrium, both cii and Aj are small, and 
we may expand cii in Aj, leading to di = LijAj. It 
is conventional to refer to di as the fiuxes, and to Ai as 
the forces. 

In the simplest case of one variable, S — dA, d — LA, 
we find R = LA^ = d'^/L, where R = J d^rTZ/T. {R is 
a function of independent variables, so it depends either 
on A or on d, but not on both simultaneously.) Clearly, 
irrespective of the tip(a), we have tip(i?) = 

The same consideration for two variables is as far as we 
need to go: Starting with S = diAi + (22^2 and employ- 
ing di = LijAj, we obtain i? as a function of either 
Ai and A2, or Ai and 02, 

R ^ LnAj + L22AI + {L12 + L2i)AiA2 (1) 
= [DAl+dl + {Li2~L2i)Aid2]/L22, (2) 

where D = L11L22 — L12L21. If R is always to be 
even under time reversal, every single term of the above 
two sums has to be even, or vanish. Now, since dS — 
Aidai, tip(S') — H-1, we have tvp{Ai) = tip(ai), or 
tip(ylid2) = — tip(AiA2), where one of the two terms 
must be odd. And since the coefficient preceding the 
odd one has to vanish, we conclude: L12 = —L21 if 
tip(AiA2) = -1, and L12 = L21 if tip(Aid2) = -1; 
or together L12 — tip{AiA2)L2i. As the subscripts 1 
and 2 are arbitrary, we have already obtained the bulk 
of the Onsager relations. 

The presence of a static magnetic field H is simple to 
include. If some kinetic coefficients depend on the field. 
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say ~ H, the coefficients are themselves odd under time 
inversion. We shall denote these odd coefficients as Jy , 
with tip(J.y ) = —1, and reserve Lij for the even ones. 
Substituting L,j by Jij in Eqs(y|2|), and reconsidering 
the tips of these terms yields two results: (i) All diago- 
nal coefficients vanish, Jn = J22 = 0. (ii) Off-diagonal 
coefficients are linked by: J12 = — tip(^iA2)^2i- This 
concludes the proof of the Onsager relation. 

It is useful to adopt a more general view of the above 
algebra, and realize two points: (i) The facts that cii is the 
time derivative of a thermodynamic variable, and Ai the 
corresponding conjugate variable, were only used once, 
to deduce their respective tips, (ii) There is nothing es- 
sential to distinguish the forces from the fluxes. Both 
can be even or odd, both vanish in equilibrium, and both 
(being of the same number) represent a complete char- 
acterization of the nonequilibrium state. Hence we may 
expand ctj in Ai, or expand Ai in ctj, or in fact expand 
a mixture of both in the rest. With these two points in 
mind, it is obvious that we could have arrived at all the 
above conclusions by starting from TS = biBi, know- 
ing only tip(5i) ~ —1, tip{Bi) — +1, and both vanish in 
equilibrium. Wc have especially 

= J2j {L'lj + Jij)Bj, with Lij = Lji, Jij = - J,;. (3) 

The generalization goes further: Any linear combina- 
tions of hi and Bi, yielding hi and Bi, will again lead 
to Eq(^) if the same conditions are met. Writing bi 
and Bi as vectors, b = (L -I- J)B, we have L = L-^, 
J = — J"^, with ^ denoting the transposition operation, 
if TS = b^B, tip(b) = -1, tip(B) = +1 hold. All four 
transformations below satisfy these three conditions, 

b = 0^b, B = e-^B, tip(©) = +1, (4) 

B = 0^b, b = ©-^B, tip(0) = -1, (5) 

b = b + *B, B:=B, (6) 

B = B-h*b, b = b, (7) 

* = tip(*) = -1, (8) 

This completes a surprisingly brief discussion of transfor- 
mations, nevertheless more general than those covered in 
the literature 0Jl^. Now, starting from Eq(||), we may 
take b = b — JB, B = B to arrive at b = LB. As L 
is symmetric, it can always be diagonalized. Dropping 
the tilde, a rescaling then reduces the dynamics to the 
strikingly simple form, b = B. 

Turning our attention now to hydrodynamic fluxes, we 
first consider heat transport in a stationary isotropic liq- 
uid with constant density, where the energy density only 
depends on the entropy density, de = Tds. (This cor- 
responds to the one- variable case above.) Defining the 
energy flux as Tf , the energy conservation e = V • (^f) 
is equivalent to 

s- V-f = (f • VT)/r. (9) 



Usually, taking f = kVT, we accept the flux f and the 
source TZ = niVT)^ as even under time reversal, oppo- 
site to s, and hence as dissipative. This is certainly not 
wrong, but it shortcuts a subtle argument and obscures 
the perfect analogy to the nonhydrodynamic case above. 
The point is, the intrinsic time inversion parity for the 
flux is negative, tip(f) = tip(s) — —1, and all three 
terms in Eq(^ are odd, with the right side characterized 
by tip(f) = -tip(VT). Such as it stands, Eq(^) is the 
exact counterpart to the starting equation of our previous 
consideration, S = a A, also a reversible equation, also 
with tip(d) = —tvp{A). Irreversibility is introduced in 
the next step, by equating these terms of different tips, 
in f = kVT as well as in a = L ^. 

The generalization to more variables is now clear: In 
all hydrodynamic theories 1^^, we have 

s - V ■ f = (^ 6,B,)/r, (10) 
i 

where again tvp{bi) — —1 and tvp{Bi) ~ -1-1, and both bi 
and Bi vanish in equilibrium. This is because the force 
(such as VT) has the tip of the variable, while the flux 
(such as f ) shares the tip with the time derivative of the 
variable. Hence the force-flux pair always has opposing 
tips, and every single previous result is duplicated. 

We start the study of nonlinear Onsager coefficients 
with the one-variable case, b — J2ni-^n + Jn)B^\ where 
the subscript „ denotes the 71*'* order. Inserting the ex- 
pression for 6 in i? = bB, we find J„ = for all n; 
switching then to _B = ^nb" we conclude Z„ = for 
even n, or B{b) — ~B{—b). 

Before going on to more variables, we need to under- 
stand one crucial point. Expanding the trace of the mo- 
mentum flux TT in the force V • 

V3 = E„^y2„+i(v-^)2"+i, (11) 

the viscosity coefficients 77271+ 1 are functions of thermo- 
dynamic variables but do not depend on the force V • t>. 
Circumstances are different if we deal with relaxing vari- 
ables, d = L2n+i^^"^^, A = dS/da. The coefficients 
L2n+i are again functions of thermodynamic variables, 
now including a, or equivalently A. The latter, how- 
ever, is also a force. As there is no way to distinguish 
A-the- variable from A-the-force, this completely under- 
mines the above result, allowing only odd orders in the 
force. Even if we admit only the linear term in the above 
sum, d = Li(A)A, we may still have terms of all orders 
in A. The electric current, j = a{E)E, is a case in point: 
The existence of a rectifier alone shows that j{E) can- 
not be an odd function of E Therefore, this type 
of consideration will yield nonlinear order results only 
for those cases in which one can clearly distinguish the 
forces from the variables, as is true for all hydrodynamic 
theories. Conversely, counter examples from relaxative 
dynamics do not disprove the time reversal invariance of 
the entropy production. 
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The insight that 72,'s invariance has something to do 
with the Casimir part of the Onsager relation has proba- 
bly crossed the mind of every practitioner of irreversible 
thermodynamics. It was Thomson who first proposed 
it. Meixner dismisses it as accidental, because he 
overlooked the difference between the hydrodynamic and 
relaxing variables, and so was misled by the fact that 
some relaxing variables are known to possess quadratic 
terms. 

Going on to find second order cross coefhcients, we 
choose the linearly diagonal representation. 



in oblivion to TZ^s inversi on p roperty, would benefit from 
embracing it, is given in |12| . 



Bi = hi + J^tji^ijk + Ji3k)bjbk, 



with Liii, L222 = 0, Lijk = Likj-, and the same for J^jk- 
First, we neglect Jijk and confine the consideration to two 
variables. There are then four independent coefficients: 
iii2 = ii2i, L122, L212 = ^221, and L211- Inserting 



the equation of motion (12) into the entropy production 



R = J2i ^i^i' obtain R = bl + Lib\b2 + 6^ + i2&i&l, 

2L■}^^) must 



where Li = 2Lii2 + L211 and L2 = ^122 '^^212 
be zero, as they precede odd terms. Substituting Bi 
for hi in i?, we find only one odd term, -6^62, the van- 
ishing coefficient of which is L211 — 2Lii2, hence L112, 
-^211 = 0. Substituting instead B2 for 62 exchanges the 
indices, leading to L122 = 2L212J and we conclude that 
all four L-coefhcients vanish. By the same procedure, we 
find that the four J-coefficients vanish, too. Next, we 
include a third variable. Three independent even coeffi- 
cients are in principle possible: L123 — Li32, L213 = L231, 
L321 = L312, in addition to the analogous odd ones. 
Again, various representations of R leave no other choice 
than to have every single coefficient vanish. Finally, we 
note that in the consideration of quadratic Onsager co- 
efficients, as there are at most three different variables 
in every term of R, we need not consider a fourth vari- 
able. This concludes the proof that no quadratic Onsager 
coefficient is compatible with the invariance of R. 

Cubic Onsager coefficients abound. Confining our at- 
tention to the two-variables case, taking Lijki to be sym- 
metric with respect to permutations of the last three in- 
dices, and again working in the linearly diagonal repre- 
sentation, the same procedure yields 6 independent co- 
efficients: 3L1112 — ^2111, 3L2221 = ^1222, -^1111, -^2222, 
L1122, and ^2112- 

The minimal conclusion one should draw from this pa- 
per is that the time inversion property of the entropy 
production TZ is worth paying much more attention to 
than is hitherto customary. This remains so even if one 
doubts the starting point of this paper, the invariance of 
TZ, because much of the physics of dissipation is codified 
in 7Z, in a highly condensed form. Bearing direct con- 
sequences on the symmetry relations between different 
kinetic coefficients, it is of great interests to those en- 
gaged in various hydrodynamic descriptions of complex 
fluids, and to those involved in kinetic theories. A spe- 
cific example of how the kinetic theory, often carried out 
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